Permittivity

 ( Electric And Magnetic Fields In Matter ) 

In electromagnetism, the absolute permittivity, often simply called permittivity and denoted by the Greek letter (epsilon), is a measure of the electric polarizability of a dielectric material. A material with high permittivity polarizes more in response to an applied electric field than a material with low permittivity, thereby storing more energy in the material. In electrostatics, the permittivity plays an important role in determining the capacitance of a capacitor.

In the simplest case, the electric displacement field resulting from an applied electric field E is

$$\backslash mathbf\{D\}\; =\; \backslash varepsilon\backslash \; \backslash mathbf\{E\}\; ~.$$

More generally, the permittivity is a thermodynamic State function.

It can depend on the frequency, Nonlinear optics, and Anisotropy of the applied field. The SI unit for permittivity is farad per meter (F/m).

The permittivity is often represented by the relative permittivity which is the ratio of the absolute permittivity and the vacuum permittivity

$$\backslash kappa\; =\; \backslash varepsilon\_\backslash mathrm\{r\}\; =\; \backslash frac\{\backslash varepsilon\}\{\backslash varepsilon\_0\}\; ~.$$

This dimensionless quantity is also often and ambiguously referred to as the permittivity. Another common term encountered for both absolute and relative permittivity is the dielectric constant which has been deprecated in physics and engineering as well as in chemistry.

By definition, a perfect vacuum has a relative permittivity of exactly 1 whereas at standard temperature and pressure, air has a relative permittivity of

Relative permittivity is directly related to electric susceptibility () by

$$\backslash chi\; =\; \backslash kappa\; -\; 1$$

otherwise written as

$$\backslash varepsilon\; =\; \backslash varepsilon\_\backslash mathrm\{r\}\backslash \; \backslash varepsilon\_0\; =\; (\; 1\; +\; \backslash chi\; )\backslash \; \backslash varepsilon\_0\; ~.$$

The term "permittivity" was introduced in the 1880s by Oliver Heaviside to complement Thomson's (1872) "permeability". Formerly written as , the designation with has been in common use since the 1950s.

---

Units The SI unit of permittivity is farad per meter (F/m or F·m⁻¹)., p. 119

$$\backslash frac\{\backslash text\{F\}\}\{\backslash text\{m\}\}\; =\; \backslash frac\{\backslash text\{C\}\}\{\backslash text\{V\}\; \{\backslash cdot\}\; \backslash text\{m\}\}\; =\; \backslash frac\{\backslash text\{C\}^2\}\{\backslash text\{N\}\; \{\backslash cdot\}\; \backslash text\{m\}^2\}\; =\; \backslash frac\{\backslash text\{C\}^2\; \{\backslash cdot\}\; \backslash text\{s\}^2\}\{\backslash text\{kg\}\; \{\backslash cdot\}\; \backslash text\{m\}^3\}=\; \backslash frac\{\backslash text\{A\}^2\; \{\backslash cdot\}\; \backslash text\{s\}^4\}\{\backslash text\{kg\}\; \{\backslash cdot\}\; \backslash text\{m\}^3\}$$

---

Explanation In electromagnetism, the electric displacement field represents the distribution of electric charges in a given medium resulting from the presence of an electric field . This distribution includes charge migration and electric dipole reorientation. Its relation to permittivity in the very simple case of linear, homogeneous, isotropic materials with "instantaneous" response to changes in electric field is:

$$\backslash mathbf\{D\}\; =\; \backslash varepsilon\backslash \; \backslash mathbf\{E\}$$

where the permittivity is a scalar. If the medium is anisotropic, the permittivity is a second rank tensor.

In general, permittivity is not a constant, as it can vary with the position in the medium, the frequency of the field applied, humidity, temperature, and other parameters. In a nonlinear optics, the permittivity can depend on the strength of the electric field. Permittivity as a function of frequency can take on real or complex values.

In SI units, permittivity is measured in farads per meter (F/m or A²·s⁴·kg⁻¹·m⁻³). The displacement field is measured in units of per square meter (C/m²), while the electric field is measured in per meter (V/m). and describe the interaction between charged objects. is related to the charge densities associated with this interaction, while is related to the forces and potential differences.

---

Vacuum permittivity The vacuum permittivity (also called permittivity of free space or the electric constant) is the ratio in Vacuum. It also appears in the Coulomb force constant,

$$k\_\backslash text\{e\}\; =\; \backslash frac\{1\}\{\backslash \; 4\backslash pi\; \backslash varepsilon\_0\backslash \; \}$$

Its value is

$$\backslash varepsilon\_0\; \backslash \; \backslash stackrel\{\backslash mathrm\{def\}\}\{=\}\backslash \; \backslash frac\{1\}\{c^2\backslash mu\_0\}\; \backslash approx\; 8.854\backslash ,187\backslash ,8128(13)\backslash times\; 10^\{-12\}\backslash text\{\; F/m\; \}$$

where

• is the speed of light in free space,
• is the vacuum permeability.

The constants and were both defined in SI units to have exact numerical values until the 2019 revision of the SI. Therefore, until that date, could be also stated exactly as a fraction, $\backslash \; \backslash tfrac\{1\}\{c^2\backslash mu\_0\}\; =\; \backslash tfrac\{1\}\{35\backslash ,950\backslash ,207\backslash ,149.472\backslash ,7056\backslash pi\}\backslash text\{\; F/m\}\backslash$ even if the result was irrational (because the fraction contained ). In contrast, the ampere was a measured quantity before 2019, but since then the ampere is now exactly defined and it is that is an experimentally measured quantity (with consequent uncertainty) and therefore so is the new 2019 definition of ( remains exactly defined before and since 2019).

---

Relative permittivity The linear permittivity of a homogeneous material is usually given relative to that of free space, as a relative permittivity (also called dielectric constant, although this term is deprecated and sometimes only refers to the static, zero-frequency relative permittivity). In an anisotropic material, the relative permittivity may be a tensor, causing birefringence. The actual permittivity is then calculated by multiplying the relative permittivity by :

$$\backslash \; \backslash varepsilon\; =\; \backslash varepsilon\_\backslash mathrm\{r\}\backslash \; \backslash varepsilon\_0\; =\; (1\; +\; \backslash chi)\backslash \; \backslash varepsilon\_0\backslash \; ,$$

where (frequently written ) is the electric susceptibility of the material.

The susceptibility is defined as the constant of proportionality (which may be a tensor) relating an electric field to the induced dielectric polarization density such that

$$\backslash \; \backslash mathbf\{P\}\backslash \; =\backslash \; \backslash varepsilon\_0\backslash \; \backslash chi\backslash \; \backslash mathbf\{E\}\backslash ;\; ,$$

where is the electric permittivity of free space.

The susceptibility of a medium is related to its relative permittivity by

$$\backslash chi\; =\; \backslash varepsilon\_\backslash mathrm\{r\}\; -\; 1\; ~.$$

So in the case of a vacuum,

$$\backslash chi\; =\; 0\; ~.$$

The susceptibility is also related to the polarizability of individual particles in the medium by the Clausius-Mossotti relation.

The electric displacement is related to the polarization density by

$$\backslash mathbf\{D\}\; =\; \backslash varepsilon\_0\backslash \; \backslash mathbf\{E\}\; +\; \backslash mathbf\{P\}\; =\; \backslash varepsilon\_0\backslash \; (1+\backslash chi)\backslash \; \backslash mathbf\{E\}\; =\; \backslash varepsilon\_\backslash mathrm\{r\}\backslash \; \backslash varepsilon\_0\backslash \; \backslash mathbf\{E\}\; ~.$$

The permittivity and permeability of a medium together determine the phase velocity of electromagnetic radiation through that medium:

$$\backslash varepsilon\; \backslash mu\; =\; \backslash frac\{\; 1\; \}\{\backslash \; v^2\; \}\; ~.$$

---

Practical applications

---

Determining capacitance The capacitance of a capacitor is based on its design and architecture, meaning it will not change with charging and discharging. The formula for capacitance in a parallel plate capacitor is written as

$$C\; =\; \backslash varepsilon\; \backslash \; \backslash frac\{A\}\{d\}$$

where $A$ is the area of one plate, $d$ is the distance between the plates, and $\backslash varepsilon$ is the permittivity of the medium between the two plates. For a capacitor with relative permittivity $\backslash kappa$, it can be said that

$$C\; =\; \backslash kappa\; \backslash \; \backslash varepsilon\_0\; \backslash frac\{A\}\{d\}$$

---

Gauss's law Permittivity is connected to electric flux (and by extension electric field) through Gauss's law. Gauss's law states that for a closed Gaussian surface, ,

$$\backslash Phi\_E\; =\; \backslash frac\{Q\_\backslash text\{enc\}\}\{\backslash varepsilon\_0\}\; =\; \backslash oint\_S\; \backslash mathbf\{E\}\; \backslash cdot\; \backslash mathrm\{d\}\; \backslash mathbf\{A\}\backslash \; ,$$ where $\backslash Phi\_E$ is the net electric flux passing through the surface, $Q\_\backslash text\{enc\}$ is the charge enclosed in the Gaussian surface, $\backslash mathbf\{E\}$ is the electric field vector at a given point on the surface, and $\backslash mathrm\{d\}\; \backslash mathbf\{A\}$ is a differential area vector on the Gaussian surface.

If the Gaussian surface uniformly encloses an insulated, symmetrical charge arrangement, the formula can be simplified to

$$E\backslash \; A\backslash \; \backslash cos\; \backslash theta\; =\; \backslash frac\{\backslash ;\; Q\_\backslash text\{enc\}\; \}\{\backslash \; \backslash varepsilon\_0\backslash \; \}\backslash \; ,$$ where $\backslash \; \backslash theta\backslash$ represents the angle between the electric field lines and the normal (perpendicular) to .

If all of the electric field lines cross the surface at 90°, the formula can be further simplified to

$$\backslash \; E\; =\; \backslash frac\{\backslash ;\; Q\_\backslash text\{enc\}\; \}\{\backslash \; \backslash varepsilon\_0\backslash \; A\backslash \; \}\; ~.$$

Because the surface area of a sphere is $\backslash \; 4\; \backslash pi\; r^2\backslash \; ,$ the electric field a distance $r$ away from a uniform, spherical charge arrangement is

$$\backslash \; E\; =\; \backslash frac\{\; Q\; \}\{\backslash \; \backslash varepsilon\_0\; A\backslash \; \}\; =\; \backslash frac\{\; Q\; \}\{\backslash \; \backslash varepsilon\_0\backslash \; \backslash left(4\backslash \; \backslash pi\backslash \; r^2\backslash right)\backslash \; \}\; =\; \backslash frac\{\; Q\; \}\{\backslash \; 4\; \backslash pi\backslash \; \backslash varepsilon\_0\backslash \; r^2\backslash \; \}\; ~.$$

This formula applies to the electric field due to a point charge, outside of a conducting sphere or shell, outside of a uniformly charged insulating sphere, or between the plates of a spherical capacitor.

---

Dispersion and causality In general, a material cannot polarize instantaneously in response to an applied field, and so the more general formulation as a function of time is

$$\backslash mathbf\{P\}(t)\; =\; \backslash varepsilon\_0\; \backslash int\_\{-\backslash infty\}^t\; \backslash chi\backslash left(t\; -\; t\text{'}\backslash right)\; \backslash mathbf\{E\}\backslash left(t\text{'}\backslash right)\; \backslash ,\; \backslash mathrm\{d\}t\text{'}\; ~.$$

That is, the polarization is a convolution of the electric field at previous times with time-dependent susceptibility given by . The upper limit of this integral can be extended to infinity as well if one defines for . An instantaneous response would correspond to a Dirac delta function susceptibility .

It is convenient to take the Fourier transform with respect to time and write this relationship as a function of frequency. Because of the convolution theorem, the integral becomes a simple product,

$$\backslash \; \backslash mathbf\{P\}(\backslash omega)\; =\; \backslash varepsilon\_0\backslash \; \backslash chi(\backslash omega)\backslash \; \backslash mathbf\{E\}(\backslash omega)\; ~.$$

This frequency dependence of the susceptibility leads to frequency dependence of the permittivity. The shape of the susceptibility with respect to frequency characterizes the dispersion properties of the material.

Moreover, the fact that the polarization can only depend on the electric field at previous times (i.e. effectively for ), a consequence of causality, imposes Kramers–Kronig constraints on the susceptibility .

---

Complex permittivity As opposed to the response of a vacuum, the response of normal materials to external fields generally depends on the frequency of the field. This frequency dependence reflects the fact that a material's polarization does not change instantaneously when an electric field is applied. The response must always be causal (arising after the applied field), which can be represented by a phase difference. For this reason, permittivity is often treated as a complex function of the (angular) frequency of the applied field:

$$\backslash varepsilon\; \backslash rightarrow\; \backslash hat\{\backslash varepsilon\}(\backslash omega)$$

(since allow specification of magnitude and phase). The definition of permittivity therefore becomes

$$D\_0\backslash \; e^\{-i\; \backslash omega\; t\}\; =\; \backslash hat\{\backslash varepsilon\}(\backslash omega)\backslash \; E\_0\backslash \; e^\{-i\; \backslash omega\; t\}\backslash \; ,$$ where

• and are the amplitudes of the displacement and electric fields, respectively,
• is the imaginary unit, .

The response of a medium to static electric fields is described by the low-frequency limit of permittivity, also called the static permittivity (also ):

$$\backslash varepsilon\_\backslash mathrm\{s\}\; =\; \backslash lim\_\{\backslash omega\; \backslash rightarrow\; 0\}\; \backslash hat\{\backslash varepsilon\}(\backslash omega)\; ~.$$

At the high-frequency limit (meaning optical frequencies), the complex permittivity is commonly referred to as (or sometimes

). At the plasma frequency and below, dielectrics behave as ideal metals, with electron gas behavior. The static permittivity is a good approximation for alternating fields of low frequencies, and as the frequency increases a measurable phase difference emerges between and . The frequency at which the phase shift becomes noticeable depends on temperature and the details of the medium. For moderate field strength (), and remain proportional, and

$$\backslash hat\{\backslash varepsilon\}\; =\; \backslash frac\{D\_0\}\{E\_0\}\; =\; |\backslash varepsilon|e^\{-i\backslash delta\}\; ~.$$

Since the response of materials to alternating fields is characterized by a complex permittivity, it is natural to separate its real and imaginary parts, which is done by convention in the following way:

$$\backslash hat\{\backslash varepsilon\}(\backslash omega)\; =\; \backslash varepsilon\text{'}(\backslash omega)\; -\; i\backslash varepsilon\text{'}\text{'}(\backslash omega)\; =\; \backslash left|\; \backslash frac\{D\_0\}\{E\_0\}\; \backslash right|\; \backslash left(\; \backslash cos\; \backslash delta\; -\; i\backslash sin\; \backslash delta\; \backslash right)\; ~.$$

where

• is the real part of the permittivity;
• is the imaginary part of the permittivity;
• is the loss angle.

The choice of sign for time-dependence, , dictates the sign convention for the imaginary part of permittivity.

The complex permittivity is usually a complicated function of frequency , since it is a superimposed description of dispersion phenomena occurring at multiple frequencies. The dielectric function must have poles only for frequencies with positive imaginary parts, and therefore satisfies the Kramers–Kronig relations. However, in the narrow frequency ranges that are often studied in practice, the permittivity can be approximated as frequency-independent or by model functions.

At a given frequency, , leads to absorption loss if it is positive (in the above sign convention) and gain if it is negative. More generally, the imaginary parts of the eigenvalues of the anisotropic dielectric tensor should be considered.

In the case of solids, the complex dielectric function is intimately connected to band structure. The primary quantity that characterizes the electronic structure of any crystalline material is the probability of photon absorption, which is directly related to the imaginary part of the optical dielectric function . The optical dielectric function is given by the fundamental expression:

$$\backslash varepsilon(\backslash omega)\; =\; 1\; +\; \backslash frac\{8\backslash pi^2\; e^2\}\{m^2\}\backslash sum\_\{c,v\}\backslash int\; W\_\{c,v\}(E)\; \backslash bigl(\; \backslash varphi\; (\backslash hbar\; \backslash omega\; -\; E)\; -\; \backslash varphi(\; \backslash hbar\backslash omega\; +\; E)\; \backslash bigr)\; \backslash ,\; \backslash mathrm\{d\}x\; ~.$$

In this expression, represents the product of the Brillouin zone-averaged transition probability at the energy with the joint density of states,

(2025). 9780750307987, Taylor and Francis. . ISBN 9780750307987

is a broadening function, representing the role of scattering in smearing out the energy levels.

In general, the broadening is intermediate between Lorentzian and Gaussian;

(1994). 9789810218645, World Scientific. . ISBN 9789810218645

Manijeh Razeghi (2025). 9780387281520, Birkhauser. . ISBN 9780387281520

for an alloy it is somewhat closer to Gaussian because of strong scattering from statistical fluctuations in the local composition on a nanometer scale.

---

Tensorial permittivity According to the Drude model of magnetized plasma, a more general expression which takes into account the interaction of the carriers with an alternating electric field at millimeter and microwave frequencies in an axially magnetized semiconductor requires the expression of the permittivity as a non-diagonal tensor:

$$\backslash mathbf\{D\}(\backslash omega)\; =\; \backslash begin\{vmatrix\}$$

   \varepsilon_1 & -i \varepsilon_2 & 0 \\
 i \varepsilon_2 &    \varepsilon_1 & 0 \\
   0             &    0             & \varepsilon_z \\
     

\end{vmatrix} \; \operatorname{\mathbf{E}}(\omega)

If vanishes, then the tensor is diagonal but not proportional to the identity and the medium is said to be a uniaxial medium, which has similar properties to a uniaxial crystal.

---

Classification of materials Materials can be classified according to their complex-valued permittivity , upon comparison of its real and imaginary components (or, equivalently, conductivity, , when accounted for in the latter). A perfect conductor has infinite conductivity, , while a perfect dielectric is a material that has no conductivity at all, ; this latter case, of real-valued permittivity (or complex-valued permittivity with zero imaginary component) is also associated with the name lossless media. Generally, when $\backslash frac\{\backslash sigma\}\{\backslash omega\; \backslash epsilon\}\; \backslash ll\; 1$ we consider the material to be a low-loss dielectric (although not exactly lossless), whereas $\backslash frac\{\backslash sigma\}\{\backslash omega\; \backslash epsilon\}\; \backslash gg\; 1$ is associated with a good conductor; such materials with non-negligible conductivity yield a large amount of dielectric loss that inhibit the propagation of electromagnetic waves, thus are also said to be lossy media. Those materials that do not fall under either limit are considered to be general media.

---

Lossy media In the case of a lossy medium, i.e. when the conduction current is not negligible, the total current density flowing is:

$$J\_\backslash text\{tot\}\backslash \; =\backslash \; J\_\backslash mathrm\{c\}\; +\; J\_\backslash mathrm\{d\}\; =\; \backslash sigma\backslash \; E\backslash \; +\backslash \; i\backslash \; \backslash omega\backslash \; \backslash varepsilon\text{'}\backslash \; E\; =\; i\backslash \; \backslash omega\backslash \; \backslash hat\{\backslash varepsilon\}\backslash \; E\backslash$$

where

• is the conductivity of the medium;
• $\backslash \; \backslash varepsilon\text{'}\backslash \; =\backslash \; \backslash varepsilon\_0\backslash \; \backslash varepsilon\_\backslash mathsf\{r\}\backslash$ is the real part of the permittivity.
• $\backslash \; \backslash hat\{\backslash varepsilon\}\backslash \; =\backslash \; \backslash varepsilon\text{'}\; -\; i\backslash \; \backslash varepsilon\text{'}\text{'}\backslash$ is the complex permittivity

Note that this is using the electrical engineering convention of the complex conjugate ambiguity; the physics/chemistry convention involves the complex conjugate of these equations.

The size of the displacement current is dependent on the frequency of the applied field ; there is no displacement current in a constant field.

In this formalism, the complex permittivity is defined as:

Kenneth L. Kaiser (2025). 9780849363726, CRC Press. . ISBN 9780849363726

$$\backslash \; \backslash hat\{\backslash varepsilon\}\backslash \; =\backslash \; \backslash varepsilon\text{'}\; \backslash left(\backslash \; 1\backslash \; -\backslash \; i\backslash \; \backslash frac\{\backslash sigma\}\{\backslash \; \backslash omega\; \backslash varepsilon\text{'}\backslash \; \}\; \backslash \; \backslash right)\backslash \; =\backslash \; \backslash varepsilon\text{'}\backslash \; -\backslash \; i\backslash \; \backslash frac\{\backslash \; \backslash sigma\backslash \; \}\{\backslash \; \backslash omega\backslash \; \}$$

In general, the absorption of electromagnetic energy by dielectrics is covered by a few different mechanisms that influence the shape of the permittivity as a function of frequency:

• First are the relaxation effects associated with permanent and induced Dipole. At low frequencies the field changes slowly enough to allow dipoles to reach before the field has measurably changed. For frequencies at which dipole orientations cannot follow the applied field because of the viscosity of the medium, absorption of the field's energy leads to energy dissipation. The mechanism of dipoles relaxing is called dielectric relaxation and for ideal dipoles is described by classic Debye relaxation.
• Second are the resonance, which arise from the rotations or vibrations of atoms, , or . These processes are observed in the neighborhood of their characteristic absorption frequencies.

The above effects often combine to cause non-linear effects within capacitors. For example, dielectric absorption refers to the inability of a capacitor that has been charged for a long time to completely discharge when briefly discharged. Although an ideal capacitor would remain at zero volts after being discharged, real capacitors will develop a small voltage, a phenomenon that is also called soakage or battery action. For some dielectrics, such as many polymer films, the resulting voltage may be less than 1–2% of the original voltage. However, it can be as much as 15–25% in the case of electrolytic capacitors or .

---

Quantum-mechanical interpretation In terms of quantum mechanics, permittivity is explained by and molecule interactions.

At low frequencies, molecules in polar dielectrics are polarized by an applied electric field, which induces periodic rotations. For example, at the microwave frequency, the microwave field causes the periodic rotation of water molecules, sufficient to break . The field does work against the bonds and the energy is absorbed by the material as heat. This is why microwave ovens work very well for materials containing water. There are two maxima of the imaginary component (the absorptive index) of water, one at the microwave frequency, and the other at far ultraviolet (UV) frequency. Both of these resonances are at higher frequencies than the operating frequency of microwave ovens.

At moderate frequencies, the energy is too high to cause rotation, yet too low to affect electrons directly, and is absorbed in the form of resonant molecular vibrations. In water, this is where the absorptive index starts to drop sharply, and the minimum of the imaginary permittivity is at the frequency of blue light (optical regime).

At high frequencies (such as UV and above), molecules cannot relax, and the energy is purely absorbed by atoms, exciting electron energy levels. Thus, these frequencies are classified as ionizing radiation.

While carrying out a complete ab initio (that is, first-principles) modelling is now computationally possible, it has not been widely applied yet. Thus, a phenomenological model is accepted as being an adequate method of capturing experimental behaviors. The Debye relaxation and the Lorentz model use a first-order and second-order (respectively) lumped system parameter linear representation (such as an RC and an LRC resonant circuit).

---

Measurement The relative permittivity of a material can be found by a variety of static electrical measurements. The complex permittivity is evaluated over a wide range of frequencies by using different variants of dielectric spectroscopy, covering nearly 21 orders of magnitude from 10⁻⁶ to 10¹⁵ hertz. Also, by using and ovens, the dielectric properties of a medium can be characterized over an array of temperatures. In order to study systems for such diverse excitation fields, a number of measurement setups are used, each adequate for a special frequency range.

Various microwave measurement techniques are outlined in Chen et al.

Typical errors for the Hakki–Coleman method employing a puck of material between conducting planes are about 0.3%.

Mailadil T. Sebastian (2025). 9780080453309, Elsevier. . ISBN 9780080453309

• Low-frequency time domain measurements ( to  Hz)
• Low-frequency frequency domain measurements ( to  Hz)
• Reflective coaxial methods ( to  Hz)
• Transmission coaxial method ( to  Hz)
• Quasi-optical methods ( to  Hz)
• Terahertz time-domain spectroscopy ( to  Hz)
• Fourier-transform methods ( to  Hz)

At infrared and optical frequencies, a common technique is ellipsometry. Dual polarisation interferometry is also used to measure the complex refractive index for very thin films at optical frequencies.

For the 3D measurement of dielectric tensors at optical frequency, Dielectric tensor tomography can be used.

---

See also

• Acoustic attenuation
• Density functional theory
• Electric-field screening
• Green–Kubo relations
• Green's function (many-body theory)
• Linear response function
• Permeability (electromagnetism)
• Rotational Brownian motion

---

Further reading

• (1973). 9780444415790, Elsevier. ISBN 9780444415790
  (volume 2 publ. 1978)
• Arthur von Hippel (2025). 9780890068038 ISBN 9780890068038
• (1966). 9780890068052 ISBN 9780890068052

---

External links

• – a chapter from an online textbook

Post Comment